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Heat flow between a large "bath" and a smaller system brings them progressively closer to thermal 
equilibrium while increasing their entropy [J]. Deviations from this trend are fluctuations involving a 
small fraction of a statistical ensemble of systems interacting with the bath: in this respect, quantum 
and classical thermodynamics are in agreement [l|, 2, HIS 01- Can there be drastic differences between 
them? Here we address a distinctly quantum mechanical setting that displays such differences: 
disturbances of thermal equilibrium between two-level systems (TLS) and a bathQ] by frequent 
and brief quantum (non-demolishing) [3, @, 0, LL0] measurements of the TLS energy-states. If the 
measurements are frequent enough to induce either the Zeno or the anti-Zeno regime, namely, the 
slowdown or speedup of the TLS relaxation [HI. \lj [l3 . \H [lq |. then the resulting entropy and 
temperature of both the system and the bath are found to be completely unrelated to what is 
expected by standard thermodynamical rules that hold for memoryless baths[2, [5|. The practical 
advantage of these anomalies is the possibility of very fast control of heat and entropy, allowing 
cooling and state-purification of quantum systems much sooner than their thermal equilibration 
time. 
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To understand the origins of the predicted anomalies, 
consider a thermal bath in equilibrium with an ensem- 
ble of quantum systems. The energy of the quantum 
systems is briefly measured. How will these systems be 
affected? Classically, their equilibrium state may remain 
intact, since measurements can be chosen to be non- 
intrusive, i.e., involve no energy exchange, but merely 
provide "snapshots" of the system. Likewise, quantum 
mechanically, nearly-ideal (projective) measurements in- 
volve no energy cost when performed by macroscopic de- 
tectors on isolated systems Q. Yet finite-time coupling, 
followed by abrupt decoupling, of two quantum ensem- 
bles, which may be viewed as the detection of one en- 
semble by the other, may cause an increase of their mean 
total energy [lg|. Here we address a different scenario per- 
taining to thermodynamics: a detector briefly (nearly- 
impulsively) measures only two-level systems (TLS) that 
are initially at thermal equilibrium with a much larger 
bath. We then ask: how will the temperature and en- 
tropy of these systems evolve both during and after the 
measurement via an interplay between the detector, the 
system and the bath? 

Such nearly-impulsive quantum measurements in the 
basis of the system (energy) eigenstates, chosen to be the 
"pointer basis" Q of the detector, must transfer energy, 
via the detector-system coupling, so as to momentarily 
interrupt (override) the system-bath interaction. This 
energy transfer, resulting in a change in system-bath en- 
tanglement, triggers the distinctly quantum dynamics of 
both the system and the bath, which subsequently re- 
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distributes their mean energy and entropy in anomalous, 
unfamiliar, ways. 

It is possible to use the information gained by such 
measurements to sort out system subensembles accord- 
ing to their measured energy, in order to extract work 
or entropy change, in the spirit of Maxwell's demon [171] . 
Here, however, we let the entire TLS ensemble evolve 
regardless of the measured result, i.e., we trace out the 
detector states. Our main concern is with the question: 
if we probe this evolution by a subsequent measurement, 
how will the outcome depend on the time-separation of 
the two? The answer should elucidate the virtually un- 
known short-time evolution of quantum systems coupled 
to a bath. 

Specifically, we analyze the scenario described above 
for a TLS with energy separation huj a , that is weakly 
coupled to a thermal bath of harmonic oscillators, char- 
acterized by a correlation (memory) time of the bath 
response t c 3> l/w a , which typically marks the onset 
of equilibrium. After equilibrium has been reached, we 
perform k = 1, ... ,K quantum non-demolition (QND) 
measurements 0, [|| 0, U£$ of the TLS energy states at 
times separated by A£& = tk+i — tk- Each measure- 
ment has a brief duration t^ -C l/u> a - Our aim is to 
explore the evolution as a function of time-separations 
between consecutive measurements in the non-Markov 
domain, Atk < l/w a <C t c . Such measurements do 
not resolve the energies of the TLS states, due to the 
time-energy uncertainty. Yet they can discriminate be- 
tween states of different symmetry, e.g., different angu- 
lar momenta. In a TLS, the evolution of the mean en- 
ergy or state populations can be identified with the ef- 
fective (spin) temperature change[18|, however rapidly it 
occurs. In this uncharted domain, we show that consecu- 
tive brief measurements entail several anomalies: (i) The 
quantum-mechanical non-commutativity of the system- 



detector and system-bath interactions causes the heat- 
up of the system at the expense of the detector-system 
coupling, but not at the expense of the coupling to the 
bath, only at very short Atfc compatible with the quan- 
tum Zeno effect fQZE)[ll|. (ii) A transition from heat- 
ing to cooling of the TLS ensemble may occur as we 
vary the interval between consecutive measurements from 
Atk *C 1/wo to Atk ~ l/wo <^ t c . This marks the tran- 
sition from Atk compatible with the QZE to those com- 
patible with the anti-Zeno effect (AZE)[H [H, [H, EH . 
Remarkably, the cooling may occur even if the bath is ini- 
tially hotter, (iii) Correspondingly, oscillations of the en- 
tropy relative to that of the equilibrium state take place, 
contrary to the Markovian notion of the second law of 

thermodynamics [2, S|- 

This scenario is governed by the following total Hamil- 
tonian of the system that interacts with the bath and is 
intermittently perturbed by the coupling of the system 
to the detector (measuring apparatus): 

H(t) = H tot + H SD (t), H tot = H s + H B +H SB . (1) 

Here Hs is the Hamiltonian of the TLS, with ground 
and excited states \g), |e), respectively; Hb is that of the 
thermal-bath composed of harmonic oscillators with en- 
ergies hoj\; Hsb = SB is the system- bath interaction 
Hamiltonian [18J] (the spin-boson interaction): a product 
of the system-dipole (or spin-flip) operator S and the 
operator B describing the bath excitations and deexci- 
tations ; and Hso(t) is the time-dependent measure- 
ment Hamiltonian that couples the system to a detec- 
tor comprised of energy-degenerate ancillae (for details 
see Supplement A). In the coupling Hamiltonians {Hsb, 
Hsd) we do not invoke the rotating- wave approximation 
(RWA) 18], namely, we do not impose energy conserva- 
tion between the system and the bath or the detector, on 
the time scales considered [15J . 

The near-equilibrium state, ptot, prior to a measure- 
ment has several pertinent characteristics (Supplement 
B): (a) It displays system-bath entanglement with off- 
diagonal matrix elements (e\p to t\g) ^ 0. (b) The sys- 
tem is described by a diagonal reduced density ma- 
trix, ps — Trsptot, in the H$ eigenbasis. (c) The 
mean interaction energy (Hsb) is negative, assuming p to t 
weakly deviates from the ground state of H to t- (Hsb) — 
(Htot) — {Hs + Hb) < 0. This comes about since the 
correction to the ground-state energy of H to t due to a 
weakly-perturbing interaction Hsb is negative (to the 
leading second order). 

We next consider the disturbance of this equilibrium 
state by a nearly-impulsive (projective) quantum mea- 
surement (t — > 0) of S, in the \g), \e) basis. The mea- 
surement correlates the TLS energy eigenstates with mu- 
tually orthogonal states of an ancillary detector and the 
latter is then averaged (traced) over. This measurement 
has distinctly quantum-mechanical consequences (Sup- 
plement A): it interrupts the system-bath interaction[13J], 
using the energy supplied by Hsd(0 < t < t) (the 
system-detector coupling) without changing (Hb) ■ It 



thus eliminates the mean system-bath interaction energy, 
whose pre-measurement value was negative, as argued 
above: 

(H SB (0)) < » (H sb (t)) = 0, (H SD (t)) = -(H SB (t)). 

(2) 
We describe the detection process as a CNOT opera- 
tion that retains the energy degeneracy of the detector 
(Hb) = 0, although its Von-Neumann entropy increases 
(Supplement A). 

After the measurement (as HsB{t ^ r) = 0), time- 
energy uncertainty at At < l/w n results in the break- 
down of the RWA, i.e., (Hs + Hb) is not conserved as 
At grows. Only (H to t) is conserved, by unitarity, un- 
til the next measurement. Hence, the post-measurement 
decrease of (Hsb) with At, signifying the restoration of 
equilibrium: (H sb (t)) = — ► (H S b(t + At)) < 0, is 
at the expense of an increase of (Hs + Hb) = (Htot) — 
(Hsb), i-e., heating of the system and the bath, com- 
bined: 



Jf [(H S ) + (H B )] 



>0, - (H SB ) 
r+At at 



\-At 



<0. (3) 



The post-measurement evolution of the system alone, 
described by ps = Trsptot, is not at all obvious. Its 
Taylor expansion holds at short evolution times, At -C 

At 2 
ps (t + At) ~ p s (r) + Atps (t) + —ps(t) + ... (4) 

The 0th order term is unchanged by the measurement, 
Ps( T ) = Ps(t ^ 0). The first derivative vanishes at t = 
"r(Ai = 0) due to the definite parity of the bath density- 
operator correlated to \g) or |e) (Supplement B). This 
initial post-measurement vanishing, ps( T ) = 0, is the 
QZE condition [Hill, El El- The time evolution of p s is 
then governed by its second time derivative ps (t) , which 
can be shown (Supplement B) to have the same sign as 
a z = |e)(e| — |.g)(g|, the population difference operator of 
the TLS. Hence, the second derivative in (J3J) is positive 
shortly after the measurement, consistently with Eq. ([3]), 
if there is no initial population inversion of the system, 
i.e., for non-negative temperature. 

The evolution of ,05 at longer times (in the 
regime of weak system-bath coupling) may be approx- 
imately described (as verified by our exact numerical 
simulations [l9j. Supplement C) by the second-order non- 
Mar kovian master equation (ME) (Fig. 1 - main panel). 
The ME for ps , on account of its diagonality, can be cast 
into the following population rate equations [15j, drop- 
ping the subscript S in what follows and setting the mea- 
surement time to be t = 0: 

Pee(t) = -f)gg(t) =Rg(t)pgg- R e (t)p ee , (5) 

/oo 
dujGT(w)sinc [(a; =F w a )t] , (6) 
-OG 

G T (u>) = (n T {u) + l)G (uj) + n T (-u})Go(-uj). (7) 



Here sinc(x) — sm ( x > ; Gt{u) is the temperature- 
dependent coupling spectrum of the bath, Gq(u>) is the 
zero-temperature coupling spectrum with peak coupling 
strength at ujq and spectral width ~ l/t c and 7Vr(w) = 
/}}_, is the inverse-temperature- (/?-) dependent popu- 
lation of bath mode u). 

The entire dynamics is determined by R e ^(t), the re- 
laxation rates of the excited (ground) states. Their non- 
Markov time-dependence yields three distinct regimes: 

(i) At short times t <C l/w "C i c the sine function in 
([5]) is much broader than Gt- The relaxation rates R e 
and R g are then equal at any temperature, indicating the 
complete breakdown of the RWA discussed above: \g) — ► 
|e) and |e) — ► |<?) transitions do not require quantum 
absorption or emission by the bath, respectively. The 
rates R e . (n) then become linear in time, manifesting the 
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This short-time regime implies the universal Zeno heating 
rate: 

J t (Pee ~ Pgg) ~ 4:R t(p gg - Pee ). (10) 

(ii) At intermediate non-Markovian times, t ~ l/w a , 
when the sine function and Gt in ([6]) have compara- 
ble widths, the relaxation rates R e ( g )(t) exhibit several 
unusual phenomena that stem from time-energy uncer- 
tainty. The change in the overlap of the sine and Gt 
functions with time results in damped aperiodic oscilla- 
tions of R e (t) and R g (t), near the frequencies loq — u> a 
and loq + w , respectively. This oscillatory time depen- 
dence that conforms neither to QZE nor to the converse 
AZE of relaxation speedup [121 \13L LL4|, will henceforth 
be dubbed the oscillatory Zeno effect (OZE). Due to the 
negativity of the sine function between its consecutive 
maxima, we can have a negative relaxation rate, which is 
completely forbidden by the RWA. Since sine [(w + u a )t] 
is much further shifted from the peak of Gt{w) than 
sine [(u> — u> a )t], R g (t) is more likely to be negative than 
R e {t) (Fig. 1(a), Fig. 2(a)). Hence, p gg {t) may grow at 
the expense of p ee (t) more than allowed by the thermal- 
equilibrium detailed balance. This may cause transient 
cooling, as detailed below. 

(iii) At long times t 3> t c , the relaxation rates attain 
their Golden-Rule (Markov) values [15| 



Re( g )(t > t e ) ~ 27tG t (±w q ). 



(11) 



The populations then approach those of an equilibrium 
Gibbs state whose temperature is equal to that of the 
thermal bath (Fig. 1 - main panel). 

We now turn to entropy dynamics. One may always 
define the entropy of ps relative to its equilibrium state 



S(ps(t)\\Po) = Tr{p s (t)\np s (t)} 



Tr{ps(t)ln-p } 

(12a) 



cr(t) = --S(p s (t)\\p ). 



(12b) 



Only in the Markovian realm, cr thus defined is identified 
as the "entropy production rate" [3, y, Q, where <r(t) ^ 
is a statement of the second law of thermodynamics in 
this realm. Since ps is diagonal, it follows (Supple- 
ment D) that cr(t) is positive iff -^ \p ee {t) — (po)ee| ^ 0, 
consistently with the interpretation of the relative en- 
tropy S(psWpo) in (|12b[) as the entropic "distance" from 
equilibrium. Conversely, whenever the oscillatory p ee {t) 
drifts away from its initial or final equilibria, cr takes 
negative values (Fig. 1(b)). 
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FIG. 1: System and bath evolution as a function of time. 
Main panel: Excited-level population as a function of time 
for initially zero-temperature product state, followed by re- 
laxation to quasi-equilibrium and then subjected to a series 
of measurements (vertical dashed lines). Measurements of 
finite duration (rjb = 0.11/w a ) (black line) results in some- 
what larger heat-up than impulsive measurements (red line), 
but the dominant effect is the same for both. Observe the 
agreement between 2 nd order master equation, two-quanta 
exchange with a discrete bath (Suppl. A), and exact numer- 
ical solution for a discrete bath of 40 modes, (a) Relaxation 
rates, R g , R e as a function of time, (b) er(t) (negative of rel- 
ative entropy rate of change), (c) Excitations as a function 
of t of the 40 modes in the two-quanta model. Parameters: 
t c = 10/ Wo, ^o = w , 7 = 0.07uj a . 

In order to realistically model the repeated measure- 
ments, i.e. give them finite duration, we assume a smooth 
temporal profile of the coupling to the detector (Supple- 
ment A). The fc-th measurement then occurs at time tu 
and has a duration of Tfe. Figure l(main panel) com- 
pares the population evolution via projective (impul- 
sive) and finite-duration measurements with t^ ~ 0.1/w a . 
Finite-duration measurements increase the Zeno-heating 
as compared to impulsive ones due to the extra energy 
supplied by the coupling to the detector. However, the 
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FIG. 2: Maximal system heating and cooling. Main panel: 
Maximum heating (lower half) and cooling (upper half) of 
the system, for different system- (a>axis) and bath- (y-axis) 
initial temperatures: as 7^ oib, where cts(B) = l/Ps(B)foua- 
(a) R g (t) and R e (t) (Eq. <j6j) ) depicted as spectral overlaps 
of the relevant functions, (b) Example of a system experi- 
encing first Zeno heating, then oscillatory-Zeno cooling, (c) 
Maximal Zeno heating (red) and subsequent maximal cooling 
(blue) as a function of common initial temperature of system 
and bath. Note the critical temperature for oscillatory-Zeno 
cooling. Parameters: t c — 10/ io a , wo = Wo/0.7, 7 = 4.36o; a . 
These effects can be strongly magnified by choosing other suit- 
able parameters, (d) Possible experimental setup. 



basic effect is seen (Fig. 1 - main panel) to be the same 
and is governed by (Hsb) change in Eq. ([3]). Counter- 
intuitively, finite-duration measurements are able to in- 
crease the cooling, despite the extra energy coupled in by 
the apparatus. 

If we repeat this procedure often enough, the TLS will 
either increasingly heat up or cool down, upon choosing 
the time intervals Atk to coincide with either peaks or 
troughs of the p ee oscillations, respectively. The mini- 
mal value of <t can also be progressively lowered with 
each measurement (Fig. 1(b)). Since consecutive mea- 
surements affect the bath and the system differently, they 
may acquire different temperatures, which then become 
the initial conditions for subsequent QZE heating or OZE 
cooling. The results are shown in Fig. 2 for both differ- 
ent (main panel) and common (Fig. 2(c)) temperatures 
of the system and the bath. Remarkably, the system 
may heat up solely due to the QZE, although the bath 
is colder, or cool down solely due to the OZE or AZE, 
although the bath is hotter (Fig. 2 - main panel). The 
bath may undergo changes in temperature and entropy 
too (Fig. 1(c)). 

One experimental realization of these effects can 
involve atoms or molecules in a microwave cavity 
(Fig. 2(d)) with controllable finite-temperature coupling 
spectrum Gt{oj) centered at luq. Measurements can be 
effected on such a TLS ensemble with resonance fre- 
quency Lu a in the microwave domain, at time intervals 
Atk ~ l/(^o ±w„), by an optical QND probe Q at fre- 
quency Up ^$> Lo a ,ojQ. The probe pulses undergo different 
Kerr-nonlinear phase shifts A<fi e or A/p g depending on 



the different symmetries (e.g., angular momenta) of |e) 
and \g). The relative abundance of A<p e and A(f> g would 
then reflect the ratio p ee (tk)/ Pgg{tk)- Such QND probing 
may be performed with time-duration much shorter than 
w„ , i.e. Lu a Tk <C 1, without resolving the energies of |e) 
and I g). 

Since non-selective measurements increase the Von- 
Neumann entropy of the detector ancillae, their entropic 
price precludes a "perpetuum mobile", if closed-cycle 
operation is attempted. Yet, if our ancillae are laser 
pulses, they are only used once and we may progressively 
change the TLS ensemble thermodynamics by consecu- 
tive pulses, disregarding their entropic or energetic price. 
The practical advantage of the predicted anomalies is 
the possibility of very rapid control of cooling and en- 
tropy, which may be attained after several measurements 
at t ^ w^ 1 and is only limited by the measurement rate. 
By contrast, conventional cooling requires much longer 
times, t>i c , to reach thermal equilibrium. 

The present findings establish a new link between fre- 
quent quantum measurements and nonequilibrium ther- 
modynamical anomalies: heat and entropy rates of 
change with the "wrong" sign, as compared to their 
usual monotonic approach to equilibrium. These anoma- 
lies are determined by the oscillatory or negative val- 
ues of the non-Markovian quantum relaxation rates at 
short times corresponding to large energy uncertainty. 
They reveal unfamiliar aspects of post-measurement 
quantum dynamics: AZE, which has been conceived as 
a means of enhancing or accelerating the initial-state 
change [12l [la. \14. fla| . here can either restore the equilib- 
rium state or further depart from it via cooling. These 
anomalies underscore the often forgotten fact that the 
system and the bath are inseparable (entangled) [a, [2fJ , 
even under weak-coupling conditions, a fact that has pro- 
found implications on their short-time dynamics. 

These results prompt further studies of a hitherto 
unexplored non-Markovian time domain where existing 
formulations of the second law of thermodynamics are 
inapplicable [2j,|5|, and neither is the common notion that 
heat always flows from hotter to colder ensembles. This 
domain may necessitate an in-depth scrutiny of quan- 
tum thermodynamical concepts. In particular, the need 
for temporal "coarse-graining" of entropy should be ex- 
amined. 
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SUPPLEMENTARY INFORMATION 

Supplement A. Measuring the Energy of a 

Two-Level System Interacting with a Bath: 

Dynamical Description 

A.l Hamiltonians and the measurement process 

We consider the following Hamiltonians: 

H(t)=H tot + H SD (t), H tot = H s + H B +H SB . (13) 

Here H to t pertains to the coupled system and bath and 
consists of: 

H s = hu a \e)(e\, (14) 

H B = k^^LO\a\ax, (15) 

A 
H SB =SB,S = \g){e\ + \e)(g\, 

B = »53(«aoa + «a4)» ( 16 ) 



where S and B are the system and bath factors, respec- 
tively, in the system-bath interaction operator, a\(a x ) 
are the annihilation (creation) operators, and k\ is the 
matrix element of the weak coupling to bath mode A. 



The detector (ancilla) qubits have energy-degenerate 
states |0)fc, |l)fe so that we may set the detector Hamil- 
tonian to be zero 



H D = 0. 



(17) 



The time-dependent system-detector coupling (to the kth 
detector) has the form 



M 



H SD (t) = J2 H so, k = J2 hk(f)\e)(e\ (I0) fefe (0| + |l> fcfc <l| - |0) fcfc <l| - |i; 

I 



fc=i 



i-fc\ 



(18) 



where 



h k (t) 



71 

4r fc 



tanh 



t-t k 

Tk 



1 



(19) 



is a smooth temporal profile of the system coupling to the 
detector qubits during the fc-th measurement that occurs 
at time tk and has a duration of Tk- 

This form of the single-measurement Hamiltonian 
HsD,k was chosen so that (dropping the index k and 
taking the measurement interval to be [0, r] for simplic- 
ity): 



e ~if^dtH SD (t)/H _ jj c _ 



(20) 



where Uc denotes to the CNOT operation (with the fcth 
detector qubit, the target qubit, denoted by the subscript 
D): 



\g)\o) D 

\e)\0) D 

1.9)11)1, 
|e)U 



D 



\g)\o) D 

\e)\l) D 
\9)\l)n 
\e)\0) D - 



(21) 



Since we take the initial state of the detector to be |0), 
only the first two rows play a role. If the measurement 



duration r is much shorter than the other time scales, 
then only H$d is non-negligible, and the entire action of 
H(t) during this time is well approximated by the CNOT 
operator Uc- This becomes exact in the impulsive limit 

T^O. 

The measurement consists in letting the TLS inter- 
act with the detector (a degenerate TLS) via Hsd- The 
measurement outcomes are averaged over (for nonselec- 
tive measurements) , by tracing out the detector degree of 
freedom. The total effect on the system density-operator 
is: 



PS >- Tr D {U c ps ® |0W0|} = \e)(e\ps\e)(e\ + \g)(g\ Ps \g)(g\ 

(22) 
i.e., the diagonal elements are unchanged, and the off- 
diagonals are erased. Since the TLS is entangled with 
the bath, the effect of the measurement in Eqs. (f2"0"l) - 
(EH) is: 



Ptot >-> Tr D {U c ptot 
= \e)(e\p tot \e)(e 



>W0|} 

\g)(g\ptot\g)(g\ 



(23) 



^e\e)(e\+pZ,\g)(g\- 



A. 2 Non-commutativity of Hsd,Hsb and the 
measurement-induced vanishing of (Hsb) 

Here we demonstrate the validity of Eq. (2) of the 
main text, based on the post-measurement increase in 
H to t following the vanishing of (Hsb)- We show that 
this effect disappears if Hsd (t) commutes with Hsb and 
Hs (and hence with H to t)- 

For the commutative case, we have: 



(Htotir)) =Tr{psBD(T)H to t} 

= Tr ^So dt[H tot+ H SD (t)] [ ptot ( )|0) 

+2f T dt'[H tot +Hsn(t')] 



(24) 



D D\ 



Htot) (25) 

Tr {ptot(0)\0) D D (0\e +l ^dt'[H tot +H SD {t')\ Htot 



xe -if T dt{H tot +H SD (t)]\. 

= Tr{p to t(0)H tot } = (H tot (0)), 



(26) 

(27) 



where psbd is the state of the combined system, bath, 
and detector. The cyclic property of the trace was used 
in (f2"6")) , and the commutativity of Hsd and H to t in (|2"T)) . 
Compare this now to our non- commutative model 
(Eqs. CE6J - (HQD) 



e -itf*ZWt)| 0)ij = ^^^ = | 1)i3 | e)(e | + |0) D |.g)( 5 |. 

(28) 
Since Hsd in Eq. (|18p commutes with Hs, we may con- 
sider the evolution of (Hsb(t)), rather than (iJ to4 (r)). 
In the impulsive limit (r — ► 0), we can drop H tot in the 
exponent of Eq. (f2"6")l . and then use the LHS of (|2"5|) to 
obtain: 



(H S B(r)) = Tr{p tot (0) 



M 



ie/Ih; 



C"SB 



(0)[/ c |0) D }. (29) 



Finally, using the RHS of ([28]) and (16|), we get: 



\UlH SB (0)Uc\0) D = -> (# SB (t)) = 0. (30) 



This expresses the vanishing of Tr {ptot{T~)HsB} due to 
the diagonality of ptot(T~) with respect to S (Suppl. B). 
Since Hd = 0, the detector mean energy is not affected 
by the CNOT action. 

Hence, the measurement-induced interruption of the 
mean interaction energy, (Hsb(t)} = 0, and the resulting 
(Hs) + (Hb) changes in Eq.(2) of the main text have a 
quantum mechanical origin: the non-commutativity of 
H S b and H S d- 



A. 3 Two-quanta approximation 



Here we assume that the system and bath, governed 
by Eqs. (TT3"|) - (TTSJ) above, were in their respective ground 
states prior to their interaction onset at t = 0, followed 
by a measurement at time tk- As shown in Suppl. B, one 
should allow for arbitrary excitations of the system and 
bath leading to an infinite hierarchy of coupled equations 
for the populations of |e), \g) and mode excitation num- 
bers. Here, for simplicity, we curtail this hierarchy, as 
is justified at short times. The wave-function driven by 
H(t) (Eq. (fl"3"|Q acquires the following form, by allowing 
the system+bath to receive or give away only or 2 exci- 
tations (the lowest two orders of the hierarchy expansion) 
through the coupling to the detector: 



\m) = E \^ l) (t))\bi) 

1=0 

H il) (t)) = Z%(t)\9)<g>\0)A + Y,^x(t)\mA®\U)A> 
A A A'jtX 

+ E4 AWl5>|2>A(8)|0>v+E^ivWI.9>|l>All>v (8) |0> A , 

A X'^X X<X' X"^X,X' 



(31) 
(32) 
(33) 



where bi is the binary representation of I, labelling the 
detector qubits. We transform to the frame where am- 



plitudes are defined by: 



a g,Q — a g,Q 


(34) 


~(l) _ -iui a t-iu x t (!) 
a e,X ~ c "e,A 


(35) 


~(l) _ -Uu x t (I) 
a g,X - e a g,X 


(36) 


„(>■) — p -i(u>x+uy )t (!) 
°g,A,A' " e a g,A,A'- 


(37) 



Using the Schrodinger equation, and integrating explic- 



itly for a { g\, 



a„ \ and a 



(0 



we obtain the following 



integro-differential matrix equation: 



the approach to equilibrium and the measurement effects, 
compared to exact numerical simulations (Suppl. C) or 
the second-order master-equation approach. 



*?>(*) 



Rx,\>(t) 



/i°(*) 



fc=i 
t 



-Ji(t)aW(t)-i/«(t)- 

M 

EM*)(aP-ai« tW) ) (38) 

1 

*A,A' E Kl„«A»e-^"(t-t')J e -- a(t'-t) (39) 

A" 

K j e fa».t+fa,*t a (0 (0) + y2 KA e^'- M ^a^ A (0) + 
£ Kve^*-*- ^41a<(0) (40) 

A^A' 



Supplement B. Bath-System Entanglement Near 
Thermal Equilibrium: Pre- and Post-Measurement 

States 



where a e (t) = {a e A } T , and Qk(l) is the decimal repre- 
sentation of {bi, ... ,1 — bk, ■■■ ,bi}, describing the flipping 
of the fcth detector qubit. This flipping occurs within the 
Hsoit) activation interval defined in Eq. (|3"8)) by hk(t) 
(Eq.©). 

Performing a brief measurement of the system at time 
tk (according to Eqs. (|38| -([40 | ) "splits" the subsequent 
evolution into two paths: (i) detection of the excited 
state with probability P^ e '(tk) = J2\ \ a e,\(tk)\ 2 , and 



» 



Wk + e) 



» 



» 



; A (t fe + e ) = a^ iA ,(t fc + e) = 0, 



where the superscript (e) denotes the excited-state out- 
come; (ii) detection of the ground state with probability 
P^)(t k ) = l~P^(t k ) anda eA (t k + e) i9) = 0, the super- 
script (g) denoting the ground-state. These subsequent 
independent evolutions destroy the system-bath correla- 
tions, and can give rise to the phenomena of Zeno heating 
and OZE cooling described in the text. 

The plots in Fig. f of the text confirm the adequacy 
of the present two-quanta approximation in describing 



The eigenstates of the total Hamiltonian H to t 
(Eq. fl!3[l) exhibit entanglement between the system and 
bath, due to the interaction term Hsb- The same is true 
of thermal states (Z^ 1 e^ ,3Htot ), at least at low tempera- 
tures and for weak coupling. 

Likewise, if the system and bath are initially in a fac- 
torizable eigenstate of H = Hs + Hb and are subse- 
quently exposed to the total Hamiltonian, including the 
interaction term, they will evolve into entangled system- 
bath states, at any temperature. If the interaction is 
turned on adiabatically, an Hq eigenstate may be ex- 
pected to evolve asymptotically into an eigenstate of 
Htot- 

The aim of this Supplement is to prove the assertion 
that in all these cases, the following properties of the 
joint state of the system and bath obtain: 

(i) ps is always diagonal (in the Hs basis), which im- 
plies in turn Tr {ptotHsB} — 0. (ii) The first moment 
of the bath excitation or deexcitation operator vanishes, 
causing the initial vanishing of the time derivative of ps 
immediately after a measurement (the Zeno effect), (iii) 
There is short-time post-measurement (Zeno) heating. 



B.l Pre-measurement evolution towards equilibrium 



We shall work in the interaction picture: 



Vi(t) = e +lHa{t - to) H S Be- lHo{t - tQ) = 



!*/(*)) 



= JH (t-t ) 



*(t)) = J7(t,*o)|*i(t )), 



^ pt ptl ptn — 1 

U(t,t Q ) = l + J2(-i) n dh dt 2 --- dt n V I (t 1 )V I (t 2 )---V I (t n ) 



EOn(t). 



(41) 
(42) 

(43) 



Let us denote the joint eigenstates of Hb and N (the where n = {rij}j and Af n is the appropriate normaliza- 
total number operator) by: tion constant. 



in) = Mi ii(4, no), 



(44) 



Consider first |^(to)) = \n) ® \g), where \n) has the 
N eigenvalue n to t- We note that V/(£)(for any i) has 
the effect of flipping the Hs state and transforming the 
bath state into a sum of states with one more or one less 
excitation. Therefore the even terms (C>2 m (t)|^(to))) in 
the perturbation expansion of |\&j(t)), Eq. (j43|) (counting 
the 1 as the zeroth term!), are superpositions of states 
with excitation numbers ntot + even, multiplied by \g), 
while the odd ones (C^m+iWI^^o))) are superpositions 
of states with nt t+odd, multiplied by |e). Let us denote 
the sum of the even terms of the series by |-B even )|e) and 
that of the odd terms as B° dd , then: 



B.2 Diagonality of ps and post-measurement 
vanishing of its first order derivative 

We wish to establish the diagonality of ps before (t ^ 
0) and after (t Js r) the measurement, and the vanishing 
of ps immediately after the measurement (t — r). Due 
to the post-measurement vanishing of the off-diagonal 
elements of ptot (Eq. (|24[) , its derivative immediately 
after the measurement, p's(r), has the form: 



PS(T) 



-i(e-^\e){g\ 



^ T \g)(e\)Tr B {B(pf g ~pg)} 



For p tot {t) = |*„ lfl (t))(*„, 9 (t)| (Eq. dUD). 



(50) 
we then 



|* 7 (t)) = U(t,t )\n,g) = \B°™(t)) ® \g) + \B° n dd (t)) ® |e) have 



= I* 



Tl.y 



(«)>• 



(45) 



Here B even (respectively, B odd ) is a sum of iV-eigenstates 
with eigenvalues differing from N by even (respectively, 
odd) numbers. 

If the initial state is of the form \^>(to)) = \n) ® |e), 
the time-evolved state is: 



(Ps) eq (t) = (e\ps(t)\g)=Tr B (e\p tot (t)\g) 



= (Bl 



(t)\B° n dd (t))=Q 



(51) 



|tt,(t)) = J7(t,to)|n,e> = |a: 

S |*n,e(*)). 



even 



(*)) ® le) 



Hence, ps is diagonal at any time t. 

In addition, we have, by virtue of (|50p at time t = t 

Ps{t) oc Tr B {6(e(. 9 )|p 4ot |e(g))} 

= (^ v f(° dd )(r)|B|B^ v r (odd) (r)) =0 (52) 



ndequation 

Now consider the initial condition that the system and 
bath are in a Gibbs state of Hq, with any inverse tem- 
perature j3: 



<e d W> ® Iff) 

(46) The same argument goes through upon permuting e «-> 
5 everywhere for p tot = |* n ,e(t)}<^n, e (<)l (Eq.®)- 

By linearity, using (fSTj) and l|5"2|) . the diagonality of 
ps(i) and the vanishing of p's immediately after the mea- 
surement are satisfied for (l48l) and (1491). 



Pt a t(t ) = Z£e- p "" =Zo'< "**Z e -, 



/3H _ 7 -l e -PHs Z -l r -l3H B 
S B 



(47) 



B.3 Second derivative positivity at arbitrary time 



z tot 7 , 



| e -/3("g+w„)| n „\/ n „i j. e _/3 ^ e+WTl '|n e)(n e | 1 For the /actoriza We thermal state, 



This initial state evolves at time t 3> to into: 

+e-«"--H-»)|* nie )(* nie |} (48) 

The resulting p to t has off-diagonal \e){g\ and |g)(e| ele- 
ments by virtue of Eqs. (|4"5)) and ([4*6")) . 

Assuming that the adiabatic theorem can be applied 
(despite the initial degeneracy of the eigenstates) , we 
have under adiabatic switching on of Hsb '■ 



Ptot(t 



^J2T, 



-0E„ 



|*n,m(*)X*r 



n m—g,e 



Ml 

E n .m being the eigenvalues of Hq. 

This state has Gibbsian form, but with the original 
Boltzmann weights. However, these weights are irrele- 
vant for the parity of the state. 



Ptot = z-i e -*ffo = z^e-^Z^e-*" 1 ', 


(53) 


we have: 




pf e = (e\p tot \e) = (e\Zj l e-^\e)Z^e-^ 


= {Ps) ee PB 


(and e <^r g). 


(54) 



For this ptot, the second derivative of ps immediately 
after the measurement is (cf. Eq. (|2~4")) ) 



(55) 



p s (r) = 2a z Tr B {B 2 (pf g -pf e )}. 

The scalar factor is positive: 



Tr B {B 2 {pf g -pf e )}=T ri 



{B 2 PB }((ps) gg -(ps) ee )>0, 
(56) 



where we have used Tr B {pf g{ee) } = {ps) gg{ee) which 
follows from the definition (Eq.(|23)l): pf e (gg) = 
(e(g)\ptot\e(g))- The first factor in (f5"6"j) is positive by 
virtue of the positivity of the operator B 2 (B being Her- 
mitian), and the second is positive iff there is no popula- 
tion inversion for the TLS. 



The combined (system- and bath-) equilibrium state 
satisfies: 



- 7-^-^tot _ 



ptot = Z e 



= 7-l p -/3O+0(fl-i B )) 



Ptot = Z e 



(57) 



Thus, for sufficiently weak coupling, Eq. (|54p dominates. 
We have made no use of the equality of the system 
and bath temperatures. The argument goes through un- 
changed for p tot = Z B 1 e- 0BHB Z- 1 e-^ Hs . 



C. Exact Numerical Simulations 

The numerical calculations have been done by 
the Multi-Configuration Time-Dependent Hartree 
(MCTDH) approach [H, [2l|. With this approach, very 
large system-bath wave functions (with more than 100 
degrees of freedom) can be propagated in time with very 
high precision. In order to treat finite temperatures we 
have sampled the Boltzmann operator in an efficient way, 
by the random-phase thermal wave function technique 



D. Entropy Dynamics 

We may always define the entropy of ps relative to 
its equilibrium state po and the negative of its rate of 
change: 



S( P s(t)\\po) 
d 



Tr{ Ps (i) lnp s (t)} ~ Tr{p s (t) lnp«$0) 

(59) 



<T(t) = --S(p s (t)\\p ). 



Only in the Markovian realm, <r thus defined is identi- 
fied as the 'entropy production rate', i.e., the net rate of 
change of the entropy of the system, after deducting the 
entropy change due to the exchange of heat Q with the 
bath at temperature T (in the absence of external work) 



3]: 



<T{t) 



dS(ps) 
dt 



1 dQ 

T~dT' 



(60) 



where S(ps) is the absolute Von-Neumann entropy of the 
system. 



It was proven by Lindblad [5|] that under any 
proper (completely positive) quantum dynamical map, 
M, the relative entropy, Eq. (|58[) cannot increase: 
S(Mp s \\Mp ) ^ S(ps\\po)- It then follows that un- 
der Markovian evolution, which is described by such a 
map, 



'(t) > 0. 



(61) 



This is a statement of the second law of 
thermodynamics [2[ . 

The quantity cr(t) defined in Eq. (|59p (which only in 
the Markovian case is identified as the entropy produc- 
tion rate of the TLS) has the following form, with (po) ee 
as the long-time (equilibrium) excitation probability 



In 



p ee (l - (p ) ee ) 

(po)ee(l ~ Pee) 



(62) 



Thus, cr(t) is positive iff ^ \p ee (t) - (po)ee| < 0, con- 
sistently with the interpretation of the relative entropy 
S(ps\\po) in (p)5|) as a kind of "distance" from equilib- 
rium. Conversely, whenever the oscillatory p ee (t) drifts 
away from its initial or final equilibria, <x must take neg- 
ative values (Fig. 1(b)). 

The condition for cr < 0, amounts to: 

R e (t)p ee -R g (t) Pgg <0 and ln^i (63) 



Lindblad's theorem on the effect of quantum maps on 
the relative entropy[5[ seems at first sight to contradict 
our results. The resolution of this apparent paradox is 
that the non-Markov evolution of ps{t) is not described 
by a map at all! This can be seen (Figs. 1,2b) from the 
fact that the density matrix oscillates and goes through 
the same value more than once, with differing subsequent 
evolution. Hence, M(0,t) : ps(0) >-> Ps(t) is a well de- 
fined completely positive map, but it is not invertible: 
different values of ps{0) may evolve into the same ps{t) 
for some particular time t. This bears likeness to the 
findings in Ref. [2(| in a different context. 
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